Magnetic miniband and magnetotransport property of a graphene superlattice 
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The eigen energy and the conductivity of a graphene sheet subject to a one-dimensional cosinu- 
soidal potential and in the presence of a magnetic field are calculated. Such a graphene superlattice 
presents three distinct magnetic miniband structures as the magnetic field increases. They are, 
respectively, the triply degenerate Landau level spectrum, the nondegenerate minibands with finite 
dispersion and the same Landau level spectrum with the pristine graphene. The ratio of the mag- 
netic length to the period of the potential function is the characteristic quantity to determine the 
electronic structure of the superlattice. Corresponding to these distinct electronic structures, the di- 
agonal conductivity presents very strong anisotropy in the weak and moderate magnetic field cases. 
But the predominant magnetotransport orientation changes from the transverse to the longitudinal 
direction of the superlattice. More interestingly, in the weak magnetic field case, the superlattice 
exhibits half-integer quantum Hall effect, but with large jump between the Ifall plateaux. Thus it 
is different from the one of the pristine graphene. 



PACS numbers: 68.65. Cd, 71.20.-b, 71.70.Di 

I. INTRODUCTION 

Since or prior to the experimental acquirement of 
graphene^, an atomically thin graphitic sheet, its 
electronic properties and potential applications were 
widely investigated. It was experimentally demon- 
strated that graphene is a gapless semiconductor 
materiali^. In the Brillouin zone of graphene, 
there are two inequivalent touching points between 
the valence and conduction bands, called the Dirac 
points^. In the low energy region relative to the 
Dirac points, the electron or hole follows a linear 
dispersion relation^i^. Thus, two massless Dirac 
fermion systems form in the vicinity of two Dirac 
points. Such an electron structure is responsible 
for most interesting electronic properties unique to 
graphene^"^. When a superlattice structure is es- 
tablished on a graphene sheet, the massless Dirac 
fermion is subject to a periodic modulation by the 
superlattice potential. Such a situation indicates 
that a graphene superlattice(GSL) perhaps exhibits 
unusual electronic characteristics, unlike both the 
pristine graphene and the ordinary semiconductor 
material^^— . Experimentally, many types of GSL 
have been fabricated and investigated. For exam- 
ple, electron-beam induced deposition of adsorbates 
on graphene membranes creates a dot array with a 
period of 5 nm^^; Epitaxially growth of graphene on 
some metallic substrates can form periodic pattern 
of supercelU^"— . These experimental progresses cer- 
tainly make the theoretical studies and device appli- 
cations of GSLs realistic. 

Accompanying the relevant experimental work, 
theoretical studies predicted many interesting elec- 
tronic properties of GSLs. For example, some re- 
cent theoretical work focused on a simple GSL 
which is constructed by exerting a one-dimensional 
strip-like or cosinusoidal periodic potential on a 
graphene sheeli^^ii^. It was found that such a 
one-dimensional GSL presents multiple Dirac points 
even in one valley of the pristine graphene^i^. The 



dependence of the number of the Dirac points on 
the strength and period of the superlattice poten- 
tial was analyzed in detail. Furthermore, when a 
vertical magnetic field is applied, the low-lying Lan- 
dau levels(LLs) become degenerate^. Accordingly, 
it was mentioned that the degeneracy of the LL ac- 
counts for the large jump of the Hall conductivity 
plateaux of the GSL^'^^. Albeit these existent theo- 
retical results, the electron structure of such a GSL 
in the presence of a magnetic field has not yet been 
comprehensively revealed. For example, magnetic 
mini-bands will be established in the GSL. But such 
a band structure and its dependence on the mag- 
netic field strength are yet unknown. The magnetic 
transport properties, namely, the diagonal and Hall 
conductivity spectra of the GSL, have not been cal- 
culated, which may be different from those of the 
pristine graphene. 

In the present work, by taking the LL states of 
the pristine graphene as the basis set, we calcu- 
late the energy spectrum of the GSL in the pres- 
ence of a magnetic field. We find that the degen- 
erate LL spectrum is only an electron structure of 
the GSL in the weak field limit, namely, in the case 
that the magnetic length is much longer than the 
period of the GSL. On the other hand, when the 
magnetic length is smaller than the period of the 
GSL, the fiat LLs evolve into magnetic minibands 
with finite dispersion. Furthermore, if the magnetic 
length is much smaller than the period of the GSL, 
i.e. the strong field case, the superlattice will re- 
store the LL spectrum of the pristine graphene. As 
for the magnetotransport properties, the diagonal 
conductivity shows notable anisotropy, especially in 
the weak magnetic field case. The Hall conductiv- 
ity shows indeed the large jumps between the ad- 
jacent plateaux in the low energy region, refiecting 
the degeneracy of the low-lying LLs. However, such 
a quantum Hall effect is destroyed in a stronger mag- 
netic field when the flat LL spectrum is replaced by 
the dispersive magnetic mini-bands. 
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The rest of the paper is organized as follows: In 
section II the theoretical method to solve the elec- 
tronic eigen-equation of the GSL in the presence of 
a vertical magnetic field is briefly elucidated. Then 
starting from the Kubo formula the conductivities 
are formulated in terms of Green function(GF). In 
section III, the numerical results about the spectra 
of the eigen energy, the density of states(DOS), the 
diagonal and Hall conductivities for the GSL are 
shown and discussed. Finally, in section IV the main 
results are briefly summarized. 



II. STRUCTURE AND THEORY 

When a magnetic field is applied perpendicular 
to the sheet of a pristine graphene, the low-energy 
electron structure in one valley, say the K valley, can 
be well described by an effective-mass Hamiltonian. 
It takes a form as- 
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where 7 is a band parameter, tt± = tTx i iTTy with 
TTa (a — x,y) being the gauged electron momen- 
tum. In Landau gauge, one has tt^ = and 
TTy — Py + eBxjc, with B denoting the magnetic 
field strength and Pa — —ihda being the momen- 
tum operator. Noting that K and K' valleys give 
the same physical result, hereafter we only treat the 
electron behavior in K valley. The Dirac point is 
taken as the zero energy. 

The eigen-equation Holnk) — e„|7ifc) can be 
solved analytically, which gives the LL spectrum of 
the pristine graphene*. It is 
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and the corresponding wavefunction of the LL state 
is given by 



and 
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where xq — Pk, is the center of the electronic har- 
monic motion; And i?n(C) is the Hermite polyno- 
mial. 

When we consider the electronic and transport 
properties of a GSL, we need to define a one- 
dimensional periodic potential along x direction, in 
addition to the Hamiltonian Hq. We choose the cos- 
inusoidal potential function given by 



V{x) = yCOs(27r2;/L). 



(7) 



It can model the periodic gate voltages applied on 
the graphene sheet. Thus, the electronic Hamilto- 
nian of such a GSL in the presence of a vertical mag- 
netic field is given by "H = Hq + V{x). Due to the 
presence of the periodic potential, it is impossible to 
obtain an analytical solution to the eigen-equation 
H\\k) — E\{k)\Xk). However, we can obtain the nu- 
merical results about it by the exact diagonalization 
technique. In doing so, we choose the LL state of 
the pristine graphene described above by Eqs(2-6) 
as the basis set. Then, we can readily find that the 
Hamiltonian matrix element is diagonal with respect 
to the electron wave vector k 

{n'k'\n\nk) = [e„5„„, + {n' k\V {x)\nk)]5kk' ■ (8) 

In the above equation, the matrix element of peri- 
odic potential can be worked out analytically. It is 
given bj^ 
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In the above two equations, hwc — It is 

d efined a s the electron cyclotron energy with I = 
^ ch/eB being the magnetic length. Lq stands for 
the linear size of the graphene sheet. And k is elec- 
tronic wave vector in the y direction. In addition, 
the normalization constant is 



1 

l/%/2 



(n = 0), 
(n^O). 



(4) 



The relevant functions are defined as 



sgn(n) = 



1 


{n > 0), 







(n = 0), 


(5) 


-1 


(71 < 0), 





. 27rxo 



exp 



±i- 



L 



Ay 
Ai! 

_^ 
2 



(10) 



with /3 = y/2nl/L, Ni = max(A,Ar') and A^2 = 
min(A^, A^'). Ln\^) stands for the associated La- 
guerre polynomials. It is defined as 

^^'(0-{^)^ie-^C^l- (11) 

The eigen wavefunction of the GSL can be for- 
mally written as 



\Xk) 
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The expansion coefficients d^/s, along with the cor- 
responding eigen energy E\{k), can be obtained by 
solving the eigen-problem of the Hamiltonian matrix 
defined above. In the actual calculation, we need to 
choose an appropriate cutofi^ index Nc which makes 
the matrix dimension finite. Of course, such a cut- 
off should guarantee the stationary results about the 
low- lying energy levels. 

By defining the retarded GF GA(i^,a;o) = 
[lo — E\{k) + irj\~^ with rj being a positive infinites- 
imal, the DOS can be expressed as 

= ^ VlmGA(w,a;o) (13) 

" "2^^^ yda;oImGA(w,xo). 

According to Kubo's linear response theory, the di- 
agonal and Hall conductivities of a two-dimensional 
electron system in the presence of a magnetic field 
can be formulated in terms of GF as 

<yo.a{e) = — 2Tr[?}aImG(£)uaImG(e)], (14) 



and 



and 



^ dujTT\v^l-mG(uj)vy — ReG(w)(15) 

ttL^ Jq duj 



—VylmG{uj)vx-j—^eG{uj)\ . 

In the above two formulae, the retarded GF operator 
is defined as G{uj) = [lj — H + iij]^^ . e denotes the 
chemical potential. Va ={i/K)[H,a] = {'y/ti)aa is 
the velocity operator which is associated with the 
Pauli matrix (Tq. When applying these formulae 
to the above GSL, by a straightforward derivation 
we can obtain the expressions of the conductivities, 
given by 
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Re[Gx{uj,xo)?lm[{Xk\vx\\'k){Xk\vy\X'k)]. 

The matrix elements of the velocity operators in- 
volved in the above two expressions are given by 



h ' nLP 



{Xk\vx\X'k) =-tlJ2diduCmCn (18) 
[sgn(r7i)(5|„j|_ij„l - sgn(7i)(5|„|j„|_i] , 



{Xk\vy\X'k) = J^d^d^G™G„ 



(19) 



[sgn(m)(5|„|_i,|„| + sgn(n)(5|„|j„|_i] . 

We can readily prove that the two relations, 
{Xk\vx\Xk) = and {Xk\vy\Xk) = {l/h)dE{k)/dk, 
always hold true for any state. These relations are 
useful for us to analyze the numerical results about 
the magnetotransport property. 



III. NUMERICAL RESULTS 

We are now in a position to perform the numer- 
ical calculations about the electronic and transport 
properties of the GSL in the presence of a magnetic 
field, according to the theoretical approach given 
above. First of all, we calculate the electronic eigen 
energy spectrum for the cases of different magnetic 
fields. We take the cutoff index Nc = 1200 which 
ensures the reliable result about the low-lying eigen 
energies close to the Dirac point. These spectra of 
a GSL are shown in Fig.l. The relevant parame- 
ters of the GSL satisfy the relation Vq — Gire^ with 
= j/L. According to previous work^'^^, such a 
GSL displays triple Dirac points in one valley of the 
pristine graphene in the case of zero magnetic field. 
We reproduce such a result by choosing the plane 
waves as the basis set to establish the Hamiltonian 
matrix of the GSL in the absence of a magnetic field. 
Such a dispersion relation is shown in Fig. 1(a). In 
Fig. 1(b) some low- lying eigen energies are shown as 
functions of the electron wavevector A:(rescaled as 
xq) for a relatively weak magnetic field. The corre- 
sponding magnetic length is larger than the period 
of the GSL. We can see that these eigen energies 
scarcely depend on xq. This is a reasonable result, 
because that the matrix elements of the Hamilto- 
nian hardly depend on the central position of the 
LL states when the magnetic length is much larger 
than the period of the GSL. Hence the eigen ener- 
gies form the flat magnetic minibands. These flat 
minibands shown in Fig. 1(b) can be viewed as the 
LL spectrum of the GSL. However, these levels are 
not exactly dispersionless, which can be clearly seen 
from the insets which blow up the weak dispersion. 
Besides, the level distribution in Fig. 1(b) is different 
from the LL spectrum of the pristine graphene. We 
can see from this figure that the three energy lev- 
els at the Dirac point are nearly degenerate. So are 
the two groups of the levels near the Dirac points. 
These results verify the LL degeneracy in the GSL 
concluded in previous work by using a distinct the- 
oretical approach2iiS. When the magnetic field in- 
creases, the flat LL spectrum changes into dispersive 
magnetic mini-bands. Such a situation is shown in 
Fig. 1(c) where the magnetic length is smaller than 
the period of the GSL by one order of magnitude. 
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FIG. 1: Low-lying eigen energies versus electronic 
wavevector k (rescaled as xo when a magnetic field is 
present) of a GSL with Vo = Quel and El ~ y/L. (a) 
The dispersion relation in the zero magnetic field case. 
The electronic wavevector along x direction is fixed at 
= 0. (b) The weak field case, i.e. I = 2.0L. The 
insets show the weak dispersion of the nearly degenerate 
levels in detail. (c) The moderate field case with I = O.IL 
and (d) the strong field case with I — O.OOIL. 



Hereafter we call such a case, namely, the occur- 
rence of the dispersive magnetic mini-bands, as the 
moderate field case. Fig. 4(d) shows the eigen energy 
spectrum when the magnetic length is far smaller 
than the period of the GSL. We can see that such a 
strong magnetic field restores the eigen energy spec- 
trum to the LL spectrum of the pristine graphene, 
described by Eq.([5]). This is because that in such a 
strong field case the periodic potential of GSL can 
be viewed as a small perturbation which only gives a 
trivial correction to the LL spectrum of the pristine 
graphene. 

In order to see the influence of the magnetic field 
on the eigen energy spectrum of the GSL more 
clearly, in Fig. 2 we show the eigen energies as func- 
tions of the magnetic length for a specific wavevec- 
tor at Xo = 0. We can see that in the moderate field 
regime the eigen energies depend sensitively on the 
variation of the magnetic field. In addition, in the 
weak field limit all the LLs tend to triply degenerate, 
which follows the triple Dirac points of the GSL. In 
Fig[3] we show the calculated DOS for such a GSL. 
The weak field case is shown in Figl31[a). The sharp 
peaks in such a DOS spectrum correspond to LLs as 
shown in Fig. 1(b). Noting that the three high peaks 
near the Dirac point arise from the triple degeneracy 
of these LLs. In the moderate field shown 
in FiglSjb), the DOS spectrum exhibits a series of 
peaks. These peaks indicate the singularities in the 
DOS, which correspond to the minima and maxima 




FIG. 2: Low- lying eigen energies versus the magnetic 
length of a GSL with Vo = Gttel at xo = 0. 
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FIG. 3: The DOS spectrum of the GSL for (a) the weak, 
(b) the moderate and (c) the strong magnetic field cases, 
respectively. 



of the magnetic mini-bands as shown in Fig. 1(c). 
In the strong field case shown in Fig. 3(c), the DOS 
spectrum is almost the same as that of the pristine 
graphene in the presence of a magnetic field. 

Now we turn to study the magnetotransport prop- 
erties of the GSL. The diagonal conductivity is cal- 
culated as a function of the chemical potential. And 
the calculated results are shown in FigUl From 
FigHIa) we can see that in the weak field case axx 
is larger than ayy. It indicates that the diagonal 
conductivity shows a strong anisotropy. We have 
known that the average velocity {Xk\vx\Xk) = for 
any eigen-state. Besides, in the weak field case the 
average velocity along y direction is ignorably small, 
namely, {Xk\vy\Xk) = {l/h)dE\{k)/dk ss 0, due to 
the flat band structure as shown in Fig. 1 (b) . There- 
fore, the diagonal conductivity in the weak field 
case is determined by the interband contributions. 
From the numerical calculation we find that the ma- 
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trix elements {Xk\vx\X'k) are notably larger than 
{Xk\vy\X'k) for any two distinct eigen-states. As a 
result, the conductivity in x direction is larger than 
that in y direction in the weak field case. In fact, the 
anisotropy of the diagonal conductivity arises from 
the anisotropic band structure around each Dirac 
points of the GSL. Previous work reported that the 
band structure of the GSL is notably anisotropic. 
Therefore, (vx) > (vy) and axx > cTyy holds true 
in zero magnetic field case^. We can infer that 
such an anisotropic band structure at each Dirac 
point of the GSL also controls the anisotropy in the 
electronic magnetotransport process presented here. 
The anisotropy in the diagonal conductivity remains 
in the moderate field case. However, opposite to the 
weak field case, the diagonal conductivity Uxx be- 
comes the smaller one. In Fig|4l[b) we can see that 
Uyy is nearly 10 times Uxx- Moreover, the peaks of 
cTxx just correspond to the dips of Uyy. We can ex- 
plain these features in the following way. The mag- 
netotransport properties of the GSL in the moderate 
field case are determined by the dispersive magnetic 
mini-bands. Then the average velocity along y direc- 
tion takes a nonzero value in any mini-band, except 
at the band-edges. Namely, {vy) — {Xk\vy\Xk) = 
{\/h)dE\{k)/dk 7^ 0. In contrast to it, the aver- 
age value of the velocity along x direction is always 
equal to zero, i.e. {Xk\vx\Xk) = in any state. Such 
a result means that Uyy has an intraband contri- 
bution, but the (Jxx does not. From Eg. p?)) . the 
formula about the diagonal conductivity, we can see 
that Uyy is mainly determined by the intraband term 
which is proportional to | (A/c|uj,|A/c) p. At the mini- 
mum and the maximum of a specific mini-band, one 
has {Xk\vy\Xk) = 0. This indicates that the dips 
in the Oyy spectrum occur at these positions. On 
the other hand, the conductivity Uxx lacks the in- 
traband contribution. Furthermore, the interband 
terms depend on the DOS of the magnetic mini- 
bands involved. This means that the conductivity 
peaks in the axx spectrum occur at the singularities 
of the DOS. And these singularities just correspond 
to the minimum or the maximum of a mini-band. 
FigHI^c) shows the diagonal conductivity spectra in 
the strong field case. We can see that the anisotropy 
disappears. And the conductivity spectrum is very 
similar to that of pristine graphene. 

The spectra of the Hall conductivity are shown 
in Figin] for the typical magnetic fields. At first, 
from Fig[5l^a) we can see that in the weak field case 
the GSL exhibits the quantum Hall effect. But the 
resolvable Hall plateaux only appear at the posi- 
tions: ±3/2(e2/ft.), ±%/2{e^/h) near the Dirac point. 
This situation is different from the half-integer quan- 
tum Hall effect of the pristine graphene. In fact, 
such a quantum Hall effect just follows the degen- 
erate LL spectrum (strictly speaking, the nearly de- 
generate LLs) in the weak field case as shown in 
Fig. 1(b). The triply degenerate low- lying LLs and 
the conduction- valence band symmetry of the energy 
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FIG. 4: (Color online) The diagonal conductivities Gxx 
and Oyy versus the chemical potential for a GSL with 
Vo = QitEl for some typical magnetic fields, (a) Mag- 
netic length I = 2.0L, (b) I = O.IL, and (c) I = O.OOIL. 
Noting that ayy in (a) is multiplied by 3 and a^x in (b) 
is multiplied by 5, in order to show these curves clearly. 



spectrum can well account for these plateaux shown 
in FiglSJ^a). When the chemical potential goes away 
from the Dirac point, the dense LLs bring about the 
oscillating character in the Hall conductivity spec- 
trum, instead of the plateau structure in the low 
energy region. In the moderate field case the fiat LL 
spectrum disappears. As a result. Hall conductiv- 
ity does not exhibit the well-defined quantum Hall 
plateaux. Such a result is shown in FiglSl^b). How- 
ever, when the magnetic field increases to the strong 
field regime, the GSL exhibits the half-integer quan- 
tum Hall effect of the pristine graphene which is 
shown in FigjSl^c). 



IV. SUMMARY 

We have studied the magnetic mini-band struc- 
ture and the magnetotransport properties of a GSL 
formed by applying a one-dimensional cosinusoidal 
potential on a graphene sheet. At first, by using the 
LL states of the pristine graphene as the basis set 
and by solving the Hamiltonian matrix of the GSL, 
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FIG. 5: The Hall conductivity axy versus the chemical 
potential for a GSL with Vb — Qtvel for some typical 
magnetic fields, (a) Magnetic length I = 2.0L, (b) I = 
O.IL, and (c) I = O.OOIL. 



mini-bands near the Dirac point. These flat mini- 
bands can be viewed as the LL spectrum of the GSL. 
Unlike the LL spectrum of the pristine graphene, 
these LLs of GSL are nearly degenerate with three 
LLs as a group. And in the weak field limit the 
LL spectrum tends to triply degenerate. When the 
magnetic length is smaller than the period of the 
GSL by one order of magnitude, namely, the mod- 
erate field case, the flat LL spectrum in the weak 
fleld case is replaced by the dispersive mini-band 
structure. And when the magnetic fleld increases 
further to enter the strong fleld regime, namely, the 
magnetic length is far shorter than period of the 
GSL, the eigen energy spectrum gets back to the LL 
spectrum of the pristine graphene. Corresponding 
to these different electronic structures, the diagonal 
and Hall conductivities present different character- 
istics in different magnetic fleld regions. In the weak 
and moderate field cases, the diagonal conductivity 
exhibits a strong anisotropy. However, the predomi- 
nant magnetotransport orientation changes from the 
transverse to the longitudinal direction of the super- 
lattice when the magnetic field increases from the 
weak to the moderate case. More interestingly, in 
the weak field case the GSL presents large jumps 
between the Hall plateaux, different from the quan- 
tum Hall effect of the pristine graphene. This feature 
arises from the degenerate LL spectrum of the GSL 
in the weak field case. 
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